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High energy emission could be produced in the interaction sites of both galactic and extragalactic 
jets with the surrounding medium. We have developed an interaction model that accounts for the 
continuous injection of relativistic electrons in the forward, reverse and recollimation shocks. We 
also performed hydrodynamical simulations to establish the physical properties in both type of 
systems. The resulting non-thermal emission is predicted assuming different values for the jet 
power, the external mass density and the source age for both FR-I galaxies and galactic micro- 
quasars. The obtained fluxes are compared to current instrument sensitivities at radio. X-ray and 
gamma-ray bands. We study in detail the connection between hard X-ray synchrotron radiation 
and gamma-ray emission and use the INTEGRAL, Fermi and AGILE, and Cherenkov telescopes 
capabilities to test it. 
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1. Introduction 

The jets of Fanaroff-Riley galaxies (of type FR-I and FR-II, [[T]]) deliver kinetic energy to the 
surrounding interstellar and intergalactic medium (ISM and IGM respectively) at a rate between 
~ 10"*^ and ~ 10"^^ erg s^^ The ejections of microquasars (/xQ), on the other hand, can also inject 
large amounts of energy into the ISM, at a level ~ 10''^ — 10^^ erg s^^ Both of jUQ and FR-I 
jet/medium interactions could be strong enough to accelerate particles and produce non-thermal 
radiation [^ ||]. Extended X-ray emission from the jets and/or lobes of the radio galaxies 3C 15 
[0], Cen A [§], Fornax A [^] and M87 [^ have been observed with Chandra. Furthermore, the 
Fermi Collaboration recently reported on the extended GeV emission from Cen A [^]. Efficient 
particle acceleration is therefore taking place at least in some FR-Is. In the case of /xQs, evidence 
of particle acceleration is found in the non-thermal radio and/or X-ray emission observed e.g. in 
SS 433 [|], XTE J1550-564 [|T§] and Cir X-1 ^), and perhaps also in GRS 1915-^105 [|T3]l and 



LS I -1-61 303 [0]. In this work we explore the leptonic non-thermal emission expected in both }xQ 
and FR-I interaction scenarios, using the results from numerical simulations coupled to a radiation 
model for the jet shocked material (reconfinement region and cocoon) and the ambient material 
shocked by the bow shock (the shell). Our results are used to make predictions for the broadband 
non-thermal fluxes (from radio to 7-rays) and we compare them with current and future instrument 
capabilities. 

2. Jet/medium interaction model 

Two symmetric jets emerge from the central source, expanding freely until they are recolli- 
mated when their ram pressure equals that of the surrounding cocoon. Further on, the jet is deceler- 
ated once the accumulated mass of the swept up ISM/ICM gas becomes similar to that carried by the 
jet. A forward shock is produced propagating into the external medium, whilst the shocked matter 
from the jet inflates the cocoon. In the FR-I scenario, the jet gets disrupted before reaching the ter- 
mination regions, and we take the jet recollimation shock as the cocoon particle accelerator. In the 
}lQ scenario the jet is assumed to remain undisrupted, and a strong reverse shock is formed near the 
jet head. In our model the mass density and pressure of the shocked regions are taken homogeneous, 
although they can evolve with time. We assume the presence of a randomly oriented magnetic field 
B in the downstream regions, derived taking the magnetic energy density to be ~ 10 % of the inter- 
nal energy density. In each shocked region, the fraction of kinetic power transferred to non-thermal 
particles is taken to be ~ 1 %. We consider the CMB and the radiation field energy density from the 
central engine (the companion star in the case of /xQ and the galaxy nucleous in the case of FR-Is), 
although the latter is only important at the very initial stages of evolution. We adopt a power-law 
spectral distribution (A^(£') = KE^, with p = 2andK such fliat jEN{E)dE = 0.1 x gjet) for the 
leptons injected at the reconfinement, bow and reverse shock fronts (we remark that the latter is 
only considered in the undisrupted jets of /xQs). Maximum energies are calculated equating the 
energy gain to synchrotron, relativistic Bremsstrahlung, Inverse Compton (IC) and adiabatic losses. 
Differently evolved populations injected all along the source age are considered at a given fgrc- For 
the synchrotron losses (?syn ?a 4 x 10'^ {B/lOjj. G)^^ (S/lTeV)^^ s), we use the magnetic field con- 
sidered above for each interaction region. Relativistic Bremsstrahlung is calculated accounting for 



Jet/medium interactions Pol Bordas 

the densities n in the downstream regions (frei.br ~ lO'^ (n/10^^cm^^)s). To compute IC losses 
(?ic ~ 1.6 X 10^^ (wrad/lO^^^erg cm^-')^^ (S/lTeV)^^ s), we consider the total radiation field en- 
ergy density Wrad- Adiabatic losses, E ss [v/r\E, are computed from the size r and the expansion 
velocity v of the emitters. Escape losses are also consideredd by taking the gyroradius of the most 



energetic particles equal to the accelerator size [15]. 



3. Hydrodynamical simulations 

Hydrodynamical simulations have been performed to further study the interaction of both jxQ 
and FR-I jets with their surroundings. A two-dimensional finite-difference code based on a high- 
resolution shock-capturing scheme has been used, which solves the equations of relativistic hy- 



drodynamics written in conservation form. The reader is referred to [|lg, [IVp for further details 
on the simulation code (see also [^). For the /xQ scenario, the jet is injected at a distance of 
10'^ cm from the compact object, with initial radius = 10'^ cm. Both the jet and the ambient 
medium are characterized with an adiabatic exponent Y = 5/3. The number density in the ambi- 
ent medium is wism = 0.3 cm^^. The velocity of the jet at injection is 0.6c, its number density 
«j = 1 .4 X 10^^ cm~^ , and its temperature T ~ 10^ ' K. These parameter values result in a jet power 
2jet = 3 X 10^^ erg s^^ For the numerical simulations of FR-I sources, the jet is injected in the 
numerical grid at 500pc from the active nucleus, with a radius of 60pc. The ambient medium 
is composed by neutral hydrogen. Its profile in pressure, density and temperature includes the 
contribution from a core region and from the galaxy group, which dominates at large distances. 
The initial jet velocity is ^ 0.87 c, its temperature ~ 4 x 10^ K, and a density and pressure ratio 
with respect to the ambient are ~ 10^^ and ~ 8, respectively, resulting in a jet kinetic luminosity 
ej = 10'*'*ergs-i. 

4. Non- thermal emission from /iQ 

Figure 1 show the spectral energy distribution (SED) for the shell, the cocoon and the jet re- 
confinement regions for the different set of parameters explored. The emission of only one jet is 
accounted for. Synchrotron emission is the channel through which the highest radiation output is 
obtained, with bolometric luminosities up to ~ 10^^ erg s^^ for powerful (Qjet = 10^' erg s^^) 
systems. At high and very high energies (HE and VHE, respectively), IC emission is the dominant 
process in the cocoon and reconfinement regions, reaching a few x 10^^" erg s^^ while in the shell 
relativistic Bremsstrahlung dominates at this energy range, with luminosities up to ~ 10^^ erg s^^ 
Notable differences are found in the reported SEDs when varying the source age ?mq from 10"* yr 
to 10^ yr. For older sources, the shell and cocoon are located at larger distances from the com- 
panion star, Mrad dccrcascs and the IC contribution gets slightly lower. Relativistic Bremsstrahlung 
emission is also larger for older sources. Higher values of ?iism make the jet to be braked at shorter 
distances from the central engine. The interaction regions have higher w^ad from the companion 
star, and the IC emission is again slightly enhanced. The relativistic Bremsstrahlung emission in 
the shell zone is also higher for denser mediums, since the luminosity is proportional to the target 
ion field density, rii ~ 4?iism- Finally, we note that in our model all the non-thermal luminosities 
scale roughly linearly with the jet power 2jet- 
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Figure 1: Spectral Energy Distribution of the non-thermal emission produced in the /iQ jet/ISM interaction. 
The contribution from the shell, the cocoon and the recollimation regions are shown in blue, red and green 
dashed lines, respectivey, whilts the overall emission is represented by the thick black continuous line. For 
the shell, synchrotron, IC and reativistic Bresstrahlung have been considered. For the recollimation and 
cocoon regions, synchrotron and IC are accounted for. Upper and lower panels correspond to a ISM particle 
density njsM = 0.1 cm^-' and 1 cm^-', respectively. Panels a), b), e), and f) show the a source age of 
fsrc = 10"*^ yr, whilst c), d), g) and h) correspond to fsrc — lO' yr. Cases a), c), e) and g) correspond to a jet 
power 2jet = 10^^ erg s^', whilst b), d), f) and h) are for Qjet = 10^^ erg s^'. See text for details on the rest 
of parameters used in the radiative model. 

5. Non-thermal emsission from FR-I sources 

The SEDs for the cocoon and the shell regions at ?sic = 10^, 3 x 10^ and 10^ yr, are shown in 
Fig. 0. The obtained radio and X-ray synchrotron luminosities in both regions are at the level of 
2 X 10"*' erg s^'. The synchrotron break frequency, corresponding to the electron energy at which 
fsyn(£^) ~ fsrc. and the highest synchrotron frequency, Vsyn max °^ ^^max' ^^^ shifted down for older 
sources. The former effect makes the radio luminosity to increase at the late stages of the evolution 
of both the cocoon and the shell, whereas the latter decreases the X-ray luminosity in the shell due 
to the decrease of Vsyn max with time. The slightly different conditions in the shell yield a higher 
break frequency, which implies a factor ~ 2 lower radio luminosity in this region compared to that 
of the cocoon. The IC luminosity grows as long as this process becomes more efficient compared to 
synchrotron and adiabatic cooling. The cocoon and the shell have similar HE luminosities, but the 
cocoon is few times brighter at VHE than the shell. In both regions the bolometric IC luminosities 
grow similarly with time, reaching ~ 10^^ and 10"^' erg s^' at HE and VHE, respectively. 

6. Discussion 



The obtained radio fluxes in a /iQ scenario would imply a flux density of ^ 150 mJy at 5 GHz 
for a source located at ~ 3 kpc. The emitting size would be of a few arcminutes, since the elec- 
tron cooling timescale is longer than the source lifetime and they can fill the whole cocoon/shell 
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Figure 2: SEDs of the non-thermal emission produced in the interaction of FR-I jets with their surroundings. 
Panels a), b) and c) correspond respectively to a source age fsrc — 10^,3 x 10^ and 10^ yr. Blue dashed lines 
show the contribution from the shell region and green dashed lines that of the cocoon. The overall emission 
is marked with a thick solid black line. 



structures. Considering this angular extension and taking a radio telescope beam size of 10", radio 
emission at a level of ~ 1 mJy beam^' could be expected. In X-rays, we find a bolometric flux in 
the range 1-10 keV of Fi-iokeV ~ 2 x 10^'^ erg s^^ cm^^. The electrons emitting at X-rays by 
synchrotron have very short time-scales, and the emitter size cannot be significantly larger than the 
accelerator itself. Although the X-rays produced in the shell tlirough relativistic Bremsstrahlung 
are expected to be quite diluted, the X-rays from the cocoon would come from a relatively small re- 
gion close to the reverse shock, and could be detectable by XMM-Newton and Chandra at scales of 
few arcseconds. At hard X-rays, INTEGRAL could detect the high energy tail of the particle popu- 
lation, although its moderate angular resolution would make difficult to resolve the sources. At HE, 
the flux between 100 MeV and 100 GeV is Fioo MeV<E<ioo GeV ~ 10^'"^ erg s^^ cm^^, below the 
Fermi sensitivity, while the integrated flux above 100 GeV is FE>iooGeV ~ 10^'^ erg s^' cm^^, also 
too low to be detectable by current Cherenkov telescopes. Taking into account the rough linearity 
between 2jet, «ism» ^src and d^^ with the gamma-ray fluxes obtained, sources with higher values of 
these quantities than the ones used here may render the /xQ jet termination regions detectable by 
current and next future gamma-ray facilities. 

For FR-I sources radio fluxes as high as ~ 10^^^ (J/100 Mpc)^^ erg cm^^ s^^ or ~ 10 Jy at 
5 GHz from a region of few times 10' (J/100 Mpc)^^ angular size are obtained. The properties of 
the radio emission from the interaction jet-medium structure are comparable with those observed 
for instance in 3c 15 |Q], in which fluxes of a few xlO^'^^ergcm^^ s^' (J ~ 300 Mpc) are found. At 
X-rays, the dominant emission comes from the cocoon, with fluxes ~ 10^'^ (J/lOO Mpc)^^ erg cm^^ : 
although limb brightening effects may increase the shell detectability. The lifetime of X-ray syn- 
chrotron electrons, ~ 10^' s, is much shorter than in radio, and ^ fsrc as well, so their radiation may 
come mostly from the inner regions of the cocoon or the bow-shock apex. The total non-thermal 
X-ray flux is roughly constant for the explored range of fgrc. although the shell contribution de- 
creases significantly with time. At HE and VHE energies, the obtained fluxes are similar for both 
the cocoon and the shell, although the latter shows a lower maximum photon energy. Gamma-ray 
emission increases with time mainly due to the increasing efficiency of the CMB IC channel. The 
obtained fluxes, for a source with fsrc ~ 10^ yr, are around ~ 10^^^ (J/lOOMpc) erg cm^^ s^^. At 
HE, such a source may require very long exposures to be detected. Nevertheless, Fermi has recently 
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detected some FR-I galaxies at a few hundred Mpc distances, including the extended radio lobes of 
Cen A (at a distance ~ 4 Mpc), presenting fluxes similar to those predicted here. At VHE, the fluxes 
would be detectable by the current instruments, although the extension of the source, of tens of ar- 
cminutes at 100 Mpc, and the steepness of the spectrum above ~ 100 GeV, may render its detection 



difficult. In the case of Cen A, detected by HESS ||18|], the emission seems to come only from the 
core, but this is expected given the large angular size of the lobes of this source, which would dilute 
its surface brightness, requiring very long observation times to render the source detectable. Fi- 
nally, we note that the fluxes showed above strongly depend on the non-thermal luminosity fraction 
going to non-thermal particles, for which we use a quite conservative value = 0.01. In the case of 
a source able to accelerate particles at a higher efficiency at the interaction shock fronts, then the 
expected non-thermal fluxes would be enhanced by a similar factor. 
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